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\0 : ABSTRACT 

O ' Multi-wavelength observations/spectroscopy of exoplanetary atmospheres are 

t^^ ■ the basis of the emerging exciting field of comparative exoplanetology. The 

HR 8799 planetary system is an ideal laboratory to study our current knowledge 
gap between massive field brown dwarfs and the cold 5-Gyr old Solar system 
planets. The HR 8799 planets have so far been imaged at J- to L-band, with 
is-^ I only upper limits available at M-band. We present here deep high-contrast Keck 

;_j ■ II adaptive optics M-band observations that show the imaging detection of 3 of 



O 



c^ 



the 4 currently known HR 8799 planets. Such detections were made possible due 
to the development of an innovative LOCI-based background subtraction scheme 
that is 3 times more efficient than a classical median background subtraction 
for Keck II AO data, representing a gain in telescope time of up to a factor 
of 9. These M-band detections extend the broad band photometric coverage out 
to ~ 5/im and provide access to the strong CO fundamental absorption band 
at 4.5yum. The new M-band photometry shows that the HR 8799 planets are 
located near the L/T-type dwarf transition, similar to what was found by other 
studies. We also confirm that the best atmospheric fits are consistent with low 
surface gravity, dusty and non-equilibrium CO/CH4 chemistry models. 

Subject headings: Instrumentation: high angular resolution, Methods: data analysis, 
Methods: observational, Techniques: image processing, : planetary systems, Planets 
and satellites: atmospheres. 1 



1. Introduction 

After more than a decade of searching, 
the direct exoplanet imaging quest was fi- 
nally successful in 2008 with the discovery 



of th r ee planetary systems flMarois et al. 



2008 



Kalas et al.l l2008l : iLagrange et al. 



20091 ) . One key advantage of this technique 
is the detection of the planet's thermal 
emission. Detailed multi-band photome- 
try and spectrometry can be acquired that 
can then be compared with atmospheric 
models to derive the planet's physical char- 
acteristics and study the effect of dust and 
molecular chemistry. With its multiple co- 
eval Jovian planets, the HR 8799 system 
is an ideal laboratory to study young plan- 
ets with low temperature/surface gravity 
atmospheres and bridge the gap between 
massive field brown dwarfs and the cold 
Solar system planets. 



H R 8799 is a 30Myr old flMarois et al. 



20101 : IZuckerman et al.l 1201 111 A5V s tar lo- 
cated 39.4pc away (JLeeuwen et al.l 120071 ) 
in the Pegasus constellation. It is clas- 
sified as__a^^l^Doradus and a A Bootis 
star (JGray et al.lll999l ). It also shows an 
IR excess (Vega-like star) consistent with 
a debris disk composed of a warm dust 
disk (6 to 15AU), a massive cold dust 
disk (90 up to 300AU) and a small dust 
particle halo extending up to a 1,000AU 
fiRhee et al.ll2007l : ISu et al.llioogh . 

Photometry and some spectroscopy of 
the HR 8799 planets are ava ilable in sev- 
eral ba nds from 1 to 3.8Atm (iMarois et al. 



2009: Metchev et al. 



Janson et al. 2010: Marois et al. 



Hinz et all 2010; Currie et al. 2011 



Barman et al.ll201lh. An ini t ial character- 



2010l:lMarois et al.ll2010l : ICurrie et all2011 : 



Barman et al.ll201ll ) has been achieved us- 



ing the available data and state-of-the-art 
atmospheric models developed for field 
brown dwarf analysis. It is clear that 
these planets are quite different than most 
field brown dwarfs, showing very dusty 
atmospheres while being cool (~1,000K) 
with no sign of methane molecular absorp- 
tion. L- and M-band detections/upper 
Umits suggest this lack of methane is 
due t o a CO/CH4 non-equilibrium chem- 



istry (ISaumon et al.l l2003l: iCushing et al. 
20061: iLeggett et al.ll2007l: Iffinz et al.ll2010l: 



Bowler et al.ll2010l : ICurrie et al.ll20 11: Ba rman et al 



2OIII ). Other teams have been able to 



fit the available photometry with patchy 
clouds an d equilibriurn chernistry atmo 
spher es (IMarois et al.l l2008l : ICurrie et al. 



2OIII ). Accurate M-band photometry can 



ization (IMarois et al.l l2008l : iBowler et al. 



help disentangle the various models, but 
it is hard on a ground-based telescope due 
to the bright thermal background mainly 
from th e telescope a nd the adaptive optics 
system (lLloydll2000l ). This background is 
particularly hard to remove at Keck, as it 
varies with time as the instrument image 
rotator (which is located very near focus) 
moves and as the AO deformable mirror 
modulates the rotator-induced background 
pattern as seen by the detector. As a 
result, conventional sky /background sub- 
traction routines leaves a spatially variable 
residual that limits final sensitivity. 

In this letter, we present the first M- 
band (4.670/im) detections of three of the 
currently four known HR 8799 planets, the 
longest wavelength at which these planets 
have been imaged. The observations are 
discussed in ^ The data reduction along 
with the new very efficient LOCI-based 



(LOCI: locally optimized combination of 
images) background subtraction technique 
are presented in f|3l Comparisons with 
field brown dwarfs and new atmospheric 
fits are shown in ^ Conclusions are sum- 
marized in §6] 

2. Observations 

The HR 8799 data was obtained on 
2009 Nov. 1 and 2 at the Keck II obser- 
vatory using the adaptive optics system, 
the NIRC2 near-infrared narrow field cam- 



era (JMcLean et al.l 120031 ) and the M-band 
filter (Ao = 4.670 /im, AA = 0.241 /im). 
No coronagraph was used for these obser- 
vations. For these two nights, 140 unsatu- 
rated images (58 for Nov. 1 and 82 for Nov. 
2) were acquired, each image consisting of 
200 coadds of 0.3s, for a total on-source 
integration time of 140min. The data 
was ta ken in angular diff erential imaging 
mode (JMarois et al.ll2006l ). In that mode, 
the instrument field rotator is adjusted 
to track the telescope pupil while leav- 
ing the field-of-view (FOV) to rotate with 
time. In addition, the telescope was nod- 
ded every few minutes between four posi- 
tions to facilitate the background subtrac- 
tion. In our data, the total FOV rotation 
is ~ 180 degrees for both nights, but both 
data sets were not acquired through transit 
(between HA [-1.6h,-0.4h] and [0.6h,l.lh] 
for Nov. 1 and [-lh,-0.6h] and [0.5h,2.2h] 
for Nov. 2). Although such observing strat- 
egy is fine for background limited data, it 
did limit the amount of speckle subtraction 
achievable at smaller separations where the 
data is still speckle noise limited and where 
planet e is now known to be located ( the 
HR 8799 transit time was actually used 
during these two nights for L-band imaging 



to confirm planet e). 

3. Data reduction 

The data reduction was performed as 
follows: we first subtract a dark image, 
divide by a fiat field and mask bad/hot 
pixels. We then remove known NIRC2 
narrow-field camer a distortions using the 
Yelda et al.l ( l2010l ) solution. The images 
are rotated by 0.252 degrees clockwise to 
put North up and the plate scale is se- 
lected to be 9.952 ± 0002mas/pixel to be 
consistent with the distortion correction of 



Yelda et al.l (120 lOf ). After these steps, we 
subtract the background (§ 13. ip and the 
speckle noise (§ 13.21) on the two data cubes 
of 58 (Nov. 1) and 82 (Nov. 2) images sepa- 
rately. We then average the two processed 
sets to produce a final image. 

3.1. LOCI-based background sub- 
traction 

The telescope was nodded between ex- 
posures to record the star at different de- 
tector locations. For each image of the se- 
quence, a classical background reduction 
would consist of subtracting the median of 
all the images where the star is sufficiently 
dithered. This standard technique is lim- 
ited because of the time variations of the 
thermal background whereas the median 
uses all the reference background images 
without distinction. As the background 
subtraction step is no different than a stan- 
dard ADI/LOCI speckle subtraction, we 
have de veloped a backgroun d subtraction 
routine (JMarois et al.ll2010bl ) based on the 
locally optimized combination of image s 
algorithm (LOCI, iLafreniere et al.l l2007l ). 
This new routine is more efficient than 
the classical median subtraction as it puts 



different weights on the background ref- 
erence images. Consider the background 
subtraction in the pth image Ip of the se- 
quence. An annulus (called D, Fig. [1]) 
is defined around the star and all refer- 
ence background images are selected where 
the star has been sufficiently dithered com- 
pared to /„. The algorithm searches for 




HR8799 



Dithered HR8799 



Fig. 1. — Image Ip where the background 
has to be subtracted (left) and the D re- 
gion inside which the LOCI coefficients are 
optimized using the reference background 
dither injages (right). 

the linear combination of these images that 
minimizes the Xp (i-G- the residual noise 
in D excluding bad pixels) which can be 
expressed as 



Xp(a 



pv-^p 



N M 



where M is the number of pixels in the im- 
age, A^ is the total number of images in the 
sequence, Xj is the jth pixel, /j is the ith. 
image of the sequence, and a^ the coeffi- 
cients to be determined. The inner radius 
of area D is chosen to avoid the bright star 
speckle noise that can bias the background 
minimization algorithm, while its outer ra- 
dius is chosen to avoid the dithered star im- 
ages. In our case, the inner and outer radii 
are 3 and 30 A/D corresponding to 24 to 
240pixels. We have then imposed a* = 



if the star is at a similar location in /j 
as in Ip. The LOC I algorithm then finds 
the optimal ali^^^* (JLafreniere et al.l 120071 : 



Marois et al.ll2010bl ) that minimizes the Xp- 
The optimized background image Bp at 
each pixel for Ip (not only in D) is then 



-Br, 



E 



a. 



fiest 



(2) 



The generated reference background im- 
age Bp is finally subtracted from Ip. 
As a*''"^** = for the images where the 
star is registered at the same location in Jj 
as in /p, Bp does not contain any flux from 
a possible companion. Moreover, if a com- 
panion exists in Jp, its impact on the a^'^*^*** 
LOCI coefficients is negligible since the 
companion flux is much smaller than the 
thermal background and D is large com- 
pared to the companion PSF. 

Once the background is subtracted, we 
sub-pixel register each image using an it- 
erative cross-correlation gaussian fit to the 
PSF core (images are unsaturated). 

3.2. Median speckle subtraction 

The speckle subtraction is performed 
following a basic API da t a redu ction as de- 
scribed in iMarois et al.l fl2006[ ). For each 
data cube, an initial reference PSF is as- 
sembled by taking the median of all images 
and is then subtracted. We then rotate the 
images to put North up and median com- 
bine them. We apply an unsharp mask to 
the resulting image (median in a 4 x 4A/D 
box) to remove the low spatial frequency 
noise and convolve it by a 0.5A/D width 
Gaussian to average the high frequency 
pixel-to-pixel noise. We finally average the 
Nov. 1 and 2 images (Fig. [2]). 
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Fig. 2. — Left panel: Final image (average of the Nov. 1 and 2 data) after subtracting the 
background noise with our new LOCI-based algorithm and the speckle noise with a basic 
median. Right panel: Same reduction, but using a classical median background subtraction. 
We applied an unsharp mask (median in a 4 x 4X/D box) on the two images and then con- 
volved them by a 0.5X/D width Gaussian. The two panels have the same linear intensity 
scale and FOV. North is up and East is left. Central panel: 3a radial contrast noise pro- 
files of our resulting median-combined images after subtracting the background with LOCI 
processing (black full line) or a median (red dashed line). Planets b, c and d fluxes are also 
plotted with la error bars along with the planet e flux upper limit. 



Without the LOCI background subtrac- 
tion, none of the planets would have been 
detected (right panel, signal-to-noise ratio 
(SNR) for d is less than 2). With the LOCI 
background subtraction, HR 8799b, c and 
d are detected (left panel, 3 to 8 SNR). 
Planet e non-detection is probably due to 
both sequences not being acquired through 
transit, thus limiting the amount of speckle 
noise being removed at small separations 
from the median subtraction. We tried 
to ap ply a more advanced LOCI algo- 



rithm (Lafreniere et al.ll2007l: iMarois et al. 



2008bl . l2010bl ) to improve the speckle re- 
duction, but as the FOV rotation ranges 
were small for both nights, no contrast gain 
was achieved. 



4. Data analysis 

Planet fluxes and positions were ob- 
tained by subtracting the planets prior 
to the speckle reduction using the stel- 
lar unsaturated PSF as the template. We 
also tried subtracting the companions prior 
to the LOCI-background algorithm and 
we have confirmed that no bias is intro- 
duced by this technique (final flux varia- 
tions smaller than 0.07%). The subtrac- 
tion was iterated by moving the planet 
template and changing its intensity until a 
minimal noise residual at the planet's loca- 
tion was achieved (inside a 1.5 A/D radius 
area; Tab. [T]for the resulting magnitudes). 
Photometric error bars were calculated 
in X/D width annulus. As expected from 



other wavelengths (JMarois et al.l |2008| ). 
HR 8799b flux is roughly a third of that 
of planets c and d. The planet's positions 
are included in Tab. [H The low SNRs and 
the large M-band PSF core result in large 
astrometric errors. A future astrometric 
HR 8799 paper using shorter wavelength 
astrometry is in preparation. 

Contrast plots (central panel of Fig. [2]) 
were obtained by calculating the noise in 
an annulus having a A/D width normal- 
ized by the stellar PSF flux (after perform- 
ing the same unsharp mask and convolu- 
tion of a 0.5 A/D Gaussian). The con- 
trast plots were then normalized at each 
separation by the estimated point source 
throughput using simulated ADI median 
process planets. Fig. [2] shows that the 
LOCI background subtraction (black full 
line) is up to ~ 3 times better than a clas- 
sical background median subtraction (red 
dotted line). If a classical background 
subtraction routine is used, an integration 
time of up to 9 times longer is required 
to reach the same LOCI background sub- 
traction contrast. This new highly efficient 
LOCI-based background subtraction rou- 
tine can be used on any data where the 
background is non-negligible and evolving 
with time. 

5. Discussions 

It has been shown that the HR 8799 
planets are an L-type extension towards 
lower effective temperatures and lower sur- 
face g r avities ( Marois et al.ll2008 ;_ Bowler et al 



2010 



Currie et al. 2011: Barman et al. 



201ll ) . The planets have also been found to 
be dusty with evidence of non-equilibrium 
CO/CH4 chemistry. If the HR 8799 plan- 
ets are plotted in K'-L' vs L'-M' dia- 



gram (Fig. [3]) agains t field brown dwarf s 
and lower mass stars (JLeggett et al.ll2007l ). 
it is found that the HR 8799 planets are 
located near the L/T-type dwarf transi- 
tion. The M-band photometry is thus 
consistent with what was found previously 
by other studies. From the field brown 
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Fig. 3. — K-L vs L-M diagram for field 
brown dwarfs, lower mass stars and the 
HR 8799 planets. Tie planets are located 
near the end of the L-type sequence, close 
to the T-type transition. 

dwarf best fi t para meters of Fig. 4-6 of 
Leggett et al.l ( 120071 ) . we also find that the 
planet L- and M-band colors are consis- 
tent wit h dusty atmospheres with vertical 
mixing. iBarman et al.l ( 1201 ll ) found a sim- 
ilar result for HR 8799b while considering 
lower surface gravities. 

If we compare ou r estimated M-band 



flux to I Currie et al.l (|201ll ) equilibrium 



chemistry model predictions, we find that 
all their best patchy-cloud atmospheric fits 



Planet 



Separation w.r.t HR8799a (E,N)" 


(1.54,0.80) 


(-0.63,0.72) 


(-0.24,-0.58) 


b 




±0.019 


±0.013 


±0.014 




Contrast Am 


10.84 ±0.30 


9.82 ±0.14 


9.44 ±0.35 


> 7.85^ 


Absolute M 


13.07 ±0.30 


12.05 ±0.14 


11.67 ±0.35 


> 10.09^ 


Absolute M fCurrie et al. 2011) 


> 11.37^^ 


> 11.22^ 


> 11.15^ 




Absolute flux (mJy) 


0.91 ±0.21 


2.33 ±0.25 


3.31 ±0.65 


< 14.23^^ 


Currie's absolute flux (mJy) 


< 4.36*^ 


< 5.00^ 


< 5.34=^ 





Table 1: HR 8799 planets M-photometry. HR 8799a apparent m agnitude is 5.21 and its 
distanc e is 39.4 ± l.Opc. Ma gnitude zeropoint is 154 Jj. liCo}Q\200u) . ^ 3cr upper limits. ^ We 
use the iMarois et all 1201 d ) astronietry for planet e. 



are rejected to ^ 3a for HR 8 799b. For 



planets c and d, ICurrie et al.l (120111 ) use 
equilibrium chemical models and focus on 
clouds to interpret the available ph o tome- 
try. The earlier work of iHinz et al.l ( 20101 ) 
did compare different chemical models to 
their M-band upper limits for all three 
planets, but did not include near-IR pho- 
tometry. Our new photometry is consis- 
tent with some models proposed by the 
two teams and will help in constraining 
their fits. 

In this letter, we compare cloudy and 
non-equilibrium model fits for all of the 
photometry for the three planets. Such 
compar ison has only been d o ne for planet b 



so fa r (iBowler et al.l l2010l : iBarman et al. 
201ll ). Using the same solar abundance 



model atmos phere grids an d proc edures 



described in iBarman et al.l (120 111 ), new 



gravities and effective temperatures are 
found (Fig. H]). For planet b, we find 
si milar model a t mosp here parameters as 



m 



Barman et al.l (120111 ) , but with a slightly 
higher gravity - log{g) = 4 rather than 3.5. 
This solar abundance model has T^s = 
HOOK and, thus, requires a very small 



radius (~ 0.7 Rj) to match the observed 
bolometric luminosity; nearly a factor of 2 
smaller than p redicted by tradition al hot- 
start models (JBaraffe et al.l |2003|). For 
planets c and d, we find temperatures and 
gravities that are fairly close to t he ex- 



pecte d values from evolution models (iMarois et al, 
20081 ) . IBarman et al.l (1201 ll ) argue that in- 
creasing the metallicity to 10 times so- 
lar could bring the Tcs/log{g) derived 
from atmosphere and evolution models for 
planet b into better agreement. The pho- 
tometry from this metal rich, evolution- 
consistent, model is also shown in Fig. HI 
This model, however, predicts an M- 
band flux that is 3.5o" brighter than ob- 
served (while all other bands agree to 2a 
or better). This may indicate an even 
higher overall metallicity and/or that a 
non-solar C-to-0 ratio is required. A 
broader exploration of the possible metal 
abundances for all four planets will be left 
to a future paper. 

6. Conclusions 

In this letter, we have estimated for the 
first time the M-band fluxes of three of 



the currently four known HR 8799 plan- 
ets. These detections were made possible 
due to the use of an innovative LOCI-based 
background subtraction routine that has 
allowed for a factor of 3 gain in contrast 
(factor of 9 in integration time) compared 
to a classical background subtraction using 
a median. This new background subtrac- 
tion routine can be used to subtract the 
background noise in any infrared data. 

We have detected HR 8799b, c and d at 
M-band from 3 to 8a. From a K'-L' and 
L'-M' color diagram, we confirm that the 
three planets are located near the end of 
the L-type sequence, close to the L- and T- 
dwarf transition region. We then derived 
new atmosphere model fits for the three 
planets. For planets c and d, temperatures 
and surface gravities are close to the ex- 
pected values from evolutionary models. 
For planet b, the solar abundance model 
fits well the broad band photometry from 1 
to 5 /xm but it requires a very small plane- 
tary radius to match the bolometric lumi- 
nosity and is thus in contradiction with the 
planet evolution models. The metal rich 
evolution-consistent model over-predicts 
the M-band flux, which may indicate an 
even higher overall metallicity and/or a 
non-solar C-to-0 ratio. Higher SNR im- 
ages will help disentangle the different 
physical and chemical parameters. 
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Fig. 4. — Comparison between the best solar-abundance model (red squares) and the broad 
band photometry (black cross) for HR 8799b (top), c (middle) and d ( bottom). For planet b , 
the photometry from the 10 times more metal rich case described in iBarman et all n201 1\) 
is plotte d in ^reen triangles . The photometric measur ements other t han for the M-band 
are from lUarois et all li200^ \20ld ): ICurrie et al] boill) . We also show lCurrie et all boill) 
upper limits for the M-band photonjetry as dotted-line. 



10 



